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ABSTRACT 

Globular clusters have long been known for presenting (at times) significant deviations from spherical sym- 
metry. While rotation has been the main proposed explanation, other complicating factors such as their constant 
interaction with the strong gravitational potential of their host galaxy have made it difficult for a consensus to 
be reached. To address this question we have obtained high-resolution spectra of WLM-l, the lone globu- 
lar cluster associated with the isolated, low-mass dwarf irregular galaxy WLM. Using archival HST WFPC2 
data, we measure the radial ellipticity profile of WLM-l, finding it to be highly elliptical, with a mean value 
of 0.17 in the region 0.5-5" — which is comparable to what is found in our Galaxy for the most elliptical 
globular clusters. There is no evidence of isophote twisting, except for the innermost regions of the cluster 
(r < 0.5"). To investigate whether the observed flattening can be ascribed to rotation, we have obtained long- 
slit high-resolution VLT/UVES spectra of this cluster along and perpendicular to the axis of flattening. Using 
cross-correlation we find that the velocity profile of the cluster is consistent with zero rotation along either axis. 
Thus neither cluster rotation nor galactic tides can be responsible for the flattened morphology of WLM- 1 . We 
argue that the required velocity dispersion anisotropy between the semi-major and semi-minor axes that would 
be required to account for the observed flattening is relatively small, of order 1 km s" 1 . Even though our errors 
preclude us from conclusively establishing that such a difference indeed exists, velocity anisotropy remains at 
present the most plausible explanation for the shape of this cluster. 
Subject headings: galaxies: star clusters — stars: kinematics 



1. INTRODUCTION 

Globular clusters have long been noted to deviat e, some- 
times substantially, from perfe ct spherical symmetry ( Shaplev 
ll917llPease & Shaple\ll917h . What is the reason why many 
globular clusters are non-spherical? A conclusive answer 
to such an important question in the study of globular clus- 
ter dynamics, with potentially important ramifications in the 
study of the evolution of globular cluster stars JNorrisll 19831 
1987), has so far not been conclusively established (see 
Mevlan & Heggiel [l997l for an extensive review and refer- 
ences). 

Several possible causes of deviations of globular cluster 
shapes from sphericity have been suggested, including tidal 
stresses due to the presence of strong galactic tidal fields 
jLongaretti & Lagoutell997llCombes et alJl9 99). gravother- 
mal shocks during th e passage of a cluster through the disk 
llKontizas et al.l fl989). velocity anisotropics, and internal ro- 
tation (which may also be acquired by interaction with a mas- 
sive galaxy; Lee et al. 2004 and references therein). In 
particular, the former two phenomena are now well docu- 
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mented, with a part icularly impressive ex ample being given 
by Palomar 5 ( Odenkir chen et all 120011) . The possibility 
that rotation may be the main dri ver of cluster ellipticities 
jRvdenlll99(Tl) is supported by the iMevlan & Ma von J198(tI) 
and Merritt, Mevlan, & Mavor ( 1997) studies of u> Centauri 
(NGC 5139), which are based upon a detailed analysis of its 
structural and dynamical profiles. Yet, it still remains unclear, 
both from an empirical and from a theoretical perspective, 
whether rotation can be the main cause of flattening of globu- 
lar star clusters. 

On the observational side, we currently lack detailed ve- 
locity profiles for most globular clusters; and even though 
we do have detailed rotational information for a few clusters, 
such as 47 Tucanae (NGC 104), a; Cen M4 (NGC 6121) and 
M15 (NGC 70781 JMavor et aljf!984l Idebhardt et aTlfl995l 
lPetersonll995llGerssen et alJ2002HZoccalil2004l) . the lack of 
a wider sample, as well as of reliable measurements of the 
variation in flattening as a function of radius for clusters en- 
compassing the full range of observed ellipticities, has hin- 
dered further progress in establishing the general correspon- 
dence, or lack thereof, between flattening and cluster rotation. 
Indeed, t he latest catalogs of g lobular cluster ellipticities, pub - 
lished bv lFrenk&FaiH dl982) and bv lWhite & Sh^im% . 
are entirely based on photographic data, consider exclusively 
the outermost cluster regions, and, relying on visual data, can 
be severely affected by differential reddening ( Ivan den Berghl 
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1982). Our group is currently building a new, detailed cat- 
alog of globular cluster ellipticities, based on 2MASS (near- 
infrar ed) data, which will greatly improve on the current situ- 
ation (Navarro, Catelan, & Stephens 2006). 

Theoretically, too, the case for rotation as one of the main 
drivers of globular cluster ellipticity is far from being settled. 
In particular, iMezal ( 120021) has recently argued that globular 
clusters can rotate to very fast rates without deviating sub- 
stantia lly from sphericity. This result is in contrast with those 
from[Alimi & Perez ( 1999), who instead were unable to pro- 
duce models of fast-rotating globular cluster-like structures 
without an accompanying deviation from sphericity, in the 
sense that their modeled structures became lengthened along 
one or two axes, depending on the initial anisotropy of the 
velocity distributions. 

An ideal test, from an empirical viewpoint, of the impact 
of rotation on the shapes of globular star clusters, is to mea- 
sure the rotation velocity of a highly flattened, isolated ex- 
tragalactic globular cluster, where there is little possibility 
that tidal effects and/or gravothermal shocks may be respon- 
sible for its flattening. The closest possible approximation to 
such an isolated, flattened system is represented by the single, 
old globular cluster in the Wolf-Lundmark-Melotte (WLM) 
dwarf irregular g alaxy (DDO 221). While distant (885 kpc; 
Hodg e et al.lll999|) . the mass of the WLM galaxy is much 
smaller than even the LMC or SMC, with M v of-14.4, -18.5, 
and -17. 1, respectively (van den Bergh 2000), thus effectively 
ruling out the possibility that the non-spherical shape of the 
WLM cluster (hereafter WLM-1) can be due to tidal effects 
or gravothermal shocks. Note, in this sense, that this is not 
nece ssarily the case even in the Magellanic Clouds (iGoodwinl 
11997b . where non-spherical globular clusters of all ages can 
be found llvan den Berghll99ll and references therein). 

In this paper we begin by reanalyzing archival HST WFPC2 
data in Section|2] where we measure the ellipticity and surface 
brightness profiles of WLM-1. We then discuss our new VLT 
UVES observations and their reduction in Section [3J and at- 
tempt to measure the rotation (Section[5} and velocity disper- 
sion (Section[6j of WLM-1. We conclude with a discussion of 
the possibility of velocity anisotropy in WLM-1 as a possible 
explanation for its observed morphology, and give some final 
remarks in Section|9] 

2. HST WFPC2 OBSERVATIONS 

In order to measure accurate ellipticity and surface bright- 
ness profiles of WLM-1, we have downloaded calibrated HST 
WFPC2 F814W and F555W images from the HST archive 
(program GO 6813, PI P. W. Hodge). In these images WLM- 
1 is centered on the PC chip. The only processing required 
was to combine the two images in each filter to reject cosmic 
rays. We also combined the images from different filters to re- 
ject the remaining few cosmic rays. This final image, rotated 
to put north at the top, is shown in Figure^ 

2.1. Cluster Ellipticity 

The ellipticity and position angle as a function of semi- 
major axis distance are measured with the Image Reduc- 
tion and Analysis Facility (IRAF 2 ) routine ELLIPSE, where 
e = l—b/a, and a and b are the semi-major and semi-minor 
axis lengths respectively. The position angle is the angle 

2 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by AURA. Inc., under cooperative agreement with the 
NSF. 
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A combined F555W/F814W HST image of WLM-1 
iHodee et al. 1999). The approximate size and orientation of the UVES slits 
are indicated. We have also inscribed ellipses illustrating the mean measured 
ellipticity (0.17) and position angle (7.1°), and the narrow extraction limits. 
The size of the image is approximately 7.7", north is up and east is left. 



(-90° < a < 90°) made by the major axis measured degrees 
east of north. The results are illustrated in Figure|2] The mean 
ellipticities and position angles for our spectral extraction re- 
gions are listed in Tabled 




0.5 1 5 

Semi-Major Axis (arcseconds) 

FIG. 2. — Ellipticity and position angle of WLM-1 as a function of semi- 
major axis distance as measured on the combined F555W and F814W frame. 
Dotted lines show Icr measurement errors. Vertical dashed lines at 0.5, 1.0, 
and 1.5" illustrate our spectral extraction limits. The horizontal dashed line 
indicates the major-axis slit position angle. 

Note that since there are many resolved stars in the HST 
image, there are many irregularities in the measurement of the 
ellipticity and position angle, especially at larger radii where 
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TABLE 1 
Ellipticity 



a(" 


) 


e 


PA(°) 


0.1 - 


0.5 


0.074 ±0.042 


4.7 ±44.1 


0.5- 


1.0 


0.153 ±0.025 


7.0±9.8 


1.0- 


1.5 


0.150±0.021 


9.4 ±6.5 


1.5- 


5.0 


0.180±0.038 


7.7 ±6.6 



these stars contribute a significant fraction of the total flux 
(see discussion in Navarro et al. 2005). Also, when the el- 
lipticity gets close to zero (nearly circular), the position angle 
becomes poorly constrained and shows large fluctuations. 

The mean ellipticity of WLM-1 in the region 0.5" < r < 
5" is 0.172 ±0.039. This ranks among the highest val- 
ues am ong Galactic globul ars according to the lists tabu- 
lated bv lFrenk&Fdi3l98l and lWhite & Shawl di987h . and 
only ~ 10 Galactic globular clusters have ellipticities greater 
than this according to the recent analysis of 2MASS data by 
Navarro et al. (2005). It has previously been noted that 
LMC and SMC globular clusters of all ages also have sig- 
nificant ellipticities, and in fact tend to be more elliptical 
than either Gal actic or M31 clusters llvan den Berghlfl99ll 
lGoodwir]ll997l and references therein), suggesting an in- 
verse relationship between galaxy mass and globular cluster 
ellipticity which would appear t o be broadly con sistent with 
the case of WLM-1. However, Harris et al. (2002) have re- 
cently shown that the gian t elliptical galaxy NGC 5128, with 
My = -20.1 mag ide Vaucouleursl [r980l. also contains a siz- 
able population of flattened clusters; 

It is also interesting to note that Han & Rvden ( 1994) have 
suggested that, while Galactic and M31 globular clusters 
shapes are broadly consistent with oblate spheroids, LMC and 
SMC clusters are instead more consistent with triaxial ellip- 
soids. As well k nown, triaxiality is often accom panied by 
isophote twisting (Williams & Schwarzschild 1979). In the 
case of WLM-1, we do find some interesting evidence that 
isophote twisting may indeed be present, but this is clearly 
restricted to the innermost 0.5" of the cluster (Fig. |2j. 

2.2. Surface Brightness Profile 

We have independently measured the surface brightness 
profile on HST WFPC2 F555W (~ V) and F814W (~ /) 
images obtained from the HST archive. Using the IRAF 
ELLIPSE task, with the ellipticity and position angle fixed 
at the mean values of 0.17 and 7.1° respectively, we mea- 
sured the average flux in each elliptical annulus, which were 
then converted to surface brightness using the photometric 
zero points (ZP) from the HST manual [ZP(555) = 22.545, 
ZP(814) = 21 .639]. The data from the F555W filter are shown 
as points in Figure|3] 

Finally w e fit the surface brightness measurements with 
iKind ( 119621) profiles using a chi-squared minimization algo- 
rithm and using the RMS scatter in each annulus as an es- 
timate of the measurement errors (although this is actually 
a measure of the true range of flux values at any particular 
semi-major axis distance rather than a measure of how accu- 
rately one can measure the surface brightness). Pixels within 
~ 0.15" of the center of the cluster are saturated, and were 
therefore excluded from the fit. These excluded points are in- 
dicated by open circles in Figure[3] The resulting fits for each 
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FIG. 3. — The background-subtracted V-band surface brightness profile of 
WLM-1 (points) and the best- fit King model (line) as determined via chi- 
squared minimization ignoring the saturated measurements (open points). 
The fit parameters, core radius (r c ) and tidal radius (r f ), are marked at the 
top, as well as the background level (dotted line) and the half-light radius 
(O,). The background-subtracted star count data from Hodge et al. 1 1999) 
have also been overplotted for comparison (histogram), with units shown on 
the right-hand axis. 



filter were then averaged together to give the final values. We 
find a best-fit core radius of r c = 0.126"±0.021", and a tidal 
radius of r, = 16.3" ± 10.4", where the magnitude of the errors 
is driven by the large scatter of fluxes in each annulus (due 
to the resolved nature of the cluster population). The fitting 
also takes into account the background surface brightness; the 
best fit values are 22.37 ±0.13 mag/arcsec 2 in F555W and 
21.45 ±0.13 mag/arcsec 2 in the F814W filter. Note that we 
have not taken into account the finite size of the HST PSF, 
which has a nominal value of le ss than 0.1" in both filters. 

Assuming (m-M)o = 24.73 (Hodge et al. 1999), one arc- 
second corresponds to 4.28 parsecs at the distance of WLM- 
1, and our measurements of the core and tidal radii convert to 
r c = 0.54 ±0.09 pc and r, = 70 ±44 pc. Both of these values 
fall squarely in the range of parameters exhibited by Galac- 
tic globulars: 0.03 < r c (pc) < 21.9 and 1.6 < r,(pc) < 214 
with median values of r c = 1.06 pc and r, = 35.4 pc (lHarrisI 
1996). This implies a very compact globular cluster, with 
c ~ 2.1 ±0.3. We cannot discard the possibility of a collapsed 
core. Note that our measured half-light r adius for the WLM-1 
cluster, along with its My = -8.74 (Sec. l6.lt . are not incon- 
sistent with the trend between r/, and My for globular clusters 
in dwarf galaxies and in the oute r Galactic halo suggested by 
Ivan den Bergh & Mackevl d2Q04t their Fig. 7). 

The structural parameters of WLM- 1 have be en previously 
derived using starcounts bv lHodge et alJ ( 119991) who found a 
core radius of 1 .09" ± 0. 14" (4.6± 0.6 pc) and a tidal radius of 
31" ± 15" (130 ±60 pc). A comparison of the background- 
subtracted star count profile and background-subtracted sur- 
face brightness profile (Fig.|5J yields a very good match down 
to ~ 2", where the star counts start to become incomplete. 
At ~ 1" the two profiles diverge; the star counts turn over 
while the surface brightness profile continues to rise down to 
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TABLE 2 
Wavelength Coverage 



Region 


Ai (A) 


A 2 (A) 


Dispersion" 


Resolution b 


Blue 


3275 


4509 


0.034 


56,700 


Red-L 


4610 


5607 


0.038 


67,220 


Red-U 


5682 


6650 


0.046 


66,450 



11 A/pixel 

b 2-pixel resolution 



~ 0.1", at which point the CCD becomes saturated. We spec- 
ulate that this incompleteness in the star counts is the source 
of Hodge's large parameter values. 

Other quantities which will be of use later in this pa- 
per include the half-light radius r/,, the central V-band sur- 
face brightness /iy(0), and the mean surface brightness in- 
side the half-light radius (nv)h- We estimate these parame- 
ters from the best-fit King model to the F555W (V) band HST 
data. We find: r h = 0.75^ arcseconds (3.2^J;| pc), ^ v (0) = 
14.83^5'4g mag/arcsec 2 , and = 17.31^;|q mag/arcsec 2 . 

3. VLT UVES OBSERVATIONS 

The new data analyzed in this paper were obtained with 
the UV- Visual Echelle Spectrograph (UVES) on the 8.2 me- 
ter VLT UT2 (Kueyen) telescope. UVES is a two-arm cross- 
dispersed echelle spectrograph with a dichroic splitting the 
red and blue components of the light. The blue light is 
recorded on a single 2K x 4K EEV CCD, and the red light 
on a mosaic of an EEV CCD and an MIT/LL CCD for the 
reddest orders. 

In order to obtain the highest spatial resolution we used a 
0.5" wide slit, which allows us to obtain velocity information 
on a similar scale. While slightly undersampling the spec- 
tral information, we chose to bin the CCDs 2 x 2 to improve 
the signal-to-noise ratio in the faint outer parts of the cluster. 
This binning yields a plate scale of 0.492"/pixel in the blue 
and 0.364"/pixel in the red. The blue slit has a length of 8" 
(~ 16 pixels) and the red slit 11" (~ 30 pixels). The wave- 
length coverage and resolution of this configuration is listed 
in TableEl 

WLM-1 (00:01:49.5, -15:27:30.7, J2000) was observed for 
3060 seconds in each of two slit positions, one aligned with 
the major axis (position angle PA = 6.8°), and another aligned 
along the mino r axis (PA = 96.8°). An HST image of the 
cluster ( iHodge et alJ l999) with the slit orientations indicated 
is shown in Figure^ 

4. UVES DATA REDUCTION 

The data reduction was carried out with IRAF. Each of 
the three UVES CCDs (Blue, Red-L, Red-U; see Table 
for wavelength coverages) were reduced independently, be- 
ginning with overscan, bias subtraction, and removal of the 
scattered inter-order light. Before dividing by the flat field, 
any variations with scales larger than a few pixels were re- 
moved from the illuminated regions by fitting and dividing by 
a high-order spline3 fit using the APFLATTEN task. The re- 
sulting "flattened" flat field was then normalized to unity and 
divided into the science frames. Lastly we located and inter- 
polated over cosmic rays and bad pixels in the science frames 
using the IRAF median filtering task CRMEDIAN. 

The center of the cluster was traced in all orders , and then 
extracted using the apertures discussed in Section 14.11 The 



calibration lamp spectrum and flat field were extracted using 
identical traces and extraction apertures as the cluster to in- 
sure the most accurate wavelength calibration. The extracted 
science spectra were then divided by the normalized extracted 
flat field to remove the blaze function and any wiggles due to 
absorption or emission features in the flat field. Finally we 
removed any large-scale problems from the extracted spectra, 
such as bad columns or cosmic ray strikes, using linear inter- 
polation from 100 pixels on either side of the blemish. 

The wavelength solution was derived from a Thorium- 
Argon calibration lamp spectrum. The solutions have an RMS 
scatter of ~ 0.003 A, and use - 1 150, - 750, and - 450 lines 
in each of the Blue, Red-L, and Red-U wavelength regions 
respectively. Note that we rederived the wavelength solution 
for each extraction. We performed one last cleaning of any re- 
maining cosmic rays in the spectra by fitting by a cubic spline 
of high enough order to fit the absorption features, and itera- 
tively replacing +4cr and -6a outliers with the fit. Finally we 
trimmed the noise-dominated ends of each order, and com- 
bined the orders into a single spectrum 

4.1. Extractions 

The key to measuring a rotation velocity in data such as 
these, where a highly concentrated object is barely resolved, 
is the careful analysis of the light coming from the edges of 
the cluster. While the center of the cluster contains most of 
the light, and hence has the highest signal-to-noise ratio, the 
slit of the spectrograph is 0.5" wide, thus we can obtain no 
velocity information in the central 0.5". The seeing, which 
ranged from 0.53" to 0.56", and the coarseness of our binned 
pixels (0.492'7pixel in the blue and 0.364"/pixel in the red) 
also blur the spatial resolution. 

We must therefore concentrate on the outer regions of the 
cluster. However, the cluster is ve ry compact, with a core ra- 
dius of 0.126" ±0.021" (Sec.O, and the light profile plum- 
mets as one moves away from the center; the cluster light 
fades into the background only ~ 3" from the center. 

We have extracted five different regions of various sizes and 
radii in hopes of obtaining the best measure of the cluster ro- 
tation profile. Since the cluster is fairly elliptical (<? ~ 0.16 in 
this region) we scale the minor-axis extraction apertures to in- 
sure that we consistently measure the same physical radius of 
the cluster. The extraction regions are listed in Table [3] The 
first region includes the entire half of the cluster, and is meant 
to give a luminosity-weighted rotation estimate. However, as 
mentioned before, due to the width of the spectrograph slit, 
the innermost regions will include an equal mix of approach- 
ing and receding stars, so the second extraction region ex- 
cludes the inner 0.5". The final three extraction regions are 
very narrow, and will be used to try to estimate the rotation 
profile of the cluster. 

When possible the extractions use background subtraction, 
where the background is calculated as the median in two re- 
gions on either side of the cluster. These background regions 
are scaled according to the ellipticity for the minor axes ex- 
tractions. Note that the slit is too narrow for background sub- 
traction on the blue chip; however, as there are no significant 
sky lines in the blue, this should not affect the results. 

One should note that the WLM- 1 clu ster is relatively large. 
The tidal radius is 16" ± 10" (Sec. l2.2> . while the slits used in 
this paper to obtain the spectra are only 8" and 11" long. Thus 
the entire slit has contributions from cluster stars, and any 
attempt at background subtraction will also subtract cluster 
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TABLE 3 
Extraction Regions' 



Major Minor 





a 


Blue 


Red 


Blue 


Red 


0.00-1.50" 
0.50-1.50" 
0.00-0.50" 
0.50-1.00" 
1.00-1.50" 
background 


0.00-3.05 
1.02-3.05 
0.00-1.02 

1.02- 2.03 

2.03- 3.05 
none 


0.00-4.12 
1.37^4.12 
0.00-1.37 
1.37-2.75 
2.75^4.12 
7.20-13.19 


0.00-2.56 
0.85-2.56 
0.00-0.85 
0.85-1.71 
1.71-2.56 
none 


0.00-3.46 
1.15-3.46 
0.00-1.15 
1.15-2.31 
2.31-3.46 
6.04-11.10 



Extraction regions in pixels. 




5160 5170 5180 

Wavelength 



FIG. 4. — Major axis spectra in three different extractions, showing the 
signal-to-noise level in each. The two spectra in each panel are from the 
north and south extractions (indicated by "N" and "S") along the major axis, 
separated by an equal vertical offset so as not to overlap in this diagram. Any 
significant rotation would be apparent as a relative shift between these two 
extractions (a 6.5 km s -1 rotation would produce a 13 km s _1 shift between 
the two spectra, corresponding to a 0.22 A shift at 5170 A). The strong fea- 
tures are Mg I lines (indicated by the vertical dotted lines). 



light. 

We plot the three narrow spectral extractions along the ma- 
jor axis in Figure|4] This very small piece of the total spec- 
trum illustrates the variation in signal-to-noise ratio with each 
extraction and the abundance of absorption features available 
for measuring velocities. 

4.2. Atmospheric Absorption 

Terrestrial atmospheric absorption lines are noticeable in 
our Red-U spectra. Two bands are especially strong: H2O 
at 5880-6020 A, and 2 at 6200-6360 A. We did not observe 
atmospheric standards, and therefore we have no way to ac- 
curately correct for this absorption. However, if not removed 
these absorption lines will bias our Red-U derived rotation 
velocity toward zero since they are at exactly the same wave- 
lengths on both sides of the cluster. 

To remove the telluric absorption lines we simply set any 
such absorption to the continuum value. We start with a high- 




6270 6280 6290 6300 6310 



wavelength 

FIG. 5. — Continuum-subtracted WLM-1 major axis spectrum before and 
after telluric masking in the wavelength region of strong atmospheric O2 ab- 
sorption. 



resolution atmospheric transmission atlas, smooth it to our 
resolution, and create a "telluric mask" for our spectra. This 
mask has regions of high transmission set to one, and absorp- 
tion regions, where the transmission is less than 95%, are set 
to a very large number. Before cross-correlating the spec- 
tra, but after continuum subtraction, we divide by this telluric 
mask, and effectively set all significant regions of atmospheric 
absorption to the continuum level. We thus lose any velocity 
information about WLM- 1 which may have been contained in 
these masked regions, but the cross-correlation will no longer 
be biased by these strong atmospheric features. A region of 
telluric absorption is shown in Figure|5]before and after mask- 
ing. 

5. CLUSTER ROTATION 

The goal of these observations is to unambiguously mea- 
sure the rotation velocity of WLM- 1 . One method of getting 
an estimate of the magnitude of this rotation would be to fol- 
low the original theoretical calculations of King ( 196JJ) for 
the ellipticity of a rotating cluster, which is equivalent to as- 
suming that the cluster is a rotating flui d body. Assuming 
a distance modulus of (m-M) = 24.73 l lHodge et a0l!999l> . 
which corresponds to 4.28 pc/arcsec, we find that the rotation 
velocity of WLM-1 should be of the order 

v rot = 6.5010ns- 1 (P—-\ I2 (—) X '\ 
VI"/ VlO 3 M /pcV Vo.16/ 

where p' is the effec tive density calculated for non- 
homogeneous clusters (van Wiik 1949). Thus if this cluster 
is indeed elliptical as a result of rotation, we might expect to 
measure a difference of ~ 13 km s _1 along the major axis, 
and zero along the minor axis (see Sec. [6]for a more detailed 
estimate of the expected rotation). 

While our spectra are not of high signal-to-noise ratio, espe- 
cially in the very narrow extractions, we do have a reasonably 

3 ftp.noao.edu/catalogs/atmospheric_transmission 
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large wavelength coverage. Thus while measuring the posi- 
tion of a single line yields a relatively large error, combining 
the results for the hundreds of lines present across the spectra 
increases the precision greatly. One way to simultaneously 
compare the positions of all features in two spectra is cross- 
correlation. 

We use th e Fou rier cross-correlation technique of 
iTonrv & Davisl (11979ft . which has been implemented in the 
IRAF FXCOR task. This routine computes the height, width, 
and location of the cross-correlation peak of two spectra. For 
each axis we correlate the north/south (major) and east/west 
(minor) extractions which lie at equal radii. The results of 
these correlations for the narrow extractions along both the 
major and minor axes are shown in Figure |6] The center of 
the peak is determined by fitting the top 5 points with Gaus- 
sian functions. This yields a direct measure of the velocity 
shift between the two, which is equal to twice the cluster ro- 
tation velocity at that distance. 
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FIG. 6. — Cross-correlation functions measured across the major (solid) 
and minor (dashed) axes, for each of the three extractions: 0-0.25" (top), 
0.25-0.50" (middle), and 0.50-75" (bottom). 

We start with the reduced, extracted spectra from each chip. 
In preparation for the cross-correlation we subtract the con- 
tinuum by fitting with a cubic spline. The telluric absorp- 
tion is removed from the Red- U spectra following the pro- 
cedure outlined in Section 1421 We then cross-correlate the 
spectra from each (Blue, Red-L, Red-U) chip. Since the re- 
ductions, extractions, and wavelength calibrations are com- 
pletely separate for each chip, this should give three indepen- 
dent, although weaker, measures of the cross-correlation ve- 
locity for each extraction. The results of this piecewise cross- 
correlation are listed in Table 0]for all symmetric extractions, 
and allow one to compare the quality of the correlation for 
each chip. We then combine the spectra from all the chips, 
and cross-correlate the result, excluding the ~ 100 A gaps 
between the chips. The results of correlating the combined 
spectra are listed in Table|5] and shown in Figure[7] 

Rather than only make one cross-correlation measurement 
to obtain the velocity at each radius, we can also cross- 
correlate all possible permutations, and find the velocities 



TABLE 4 

Piecewise Cross-Correlation Velocities 3 



region 


0.0-0.5" 


0.5-1.0" 


1.0-1.5" 


\4 A TOD AYIC 

MAJ UK-AAlo 








Blue 


+0.09 ±0.16 


+1.77 ±0.55 


-0.39 ±1.58 


Red-1 


-0.20 ±0.1 5 


-0.28 ±0.78 


+0.38 ±0.57 


Red-u 


+0.02 ±0.11 


-0.23 ± 0.73 


-0.20 ±0.40 


Combined 


-0.09 ± 0.08 


-0.17 ±0.41 


-0.17 ±0.27 


MINOR-AXIS 








Blue 


-0.02 ±0.15 


-0.38 ±0.52 


-1.48 ±0.79 


Red-L 


+0.05 ±0.16 


+0.02 ± 0.65 


-0.06 ±0.54 


Red-U 


+0.03 ±0.11 


-0.96 ±0.65 


-0.45 ±0.49 


Combined 


-0.03 ±0.08 


-0.56 ±0.32 


-0.39 ±0.31 



"Velocities in km s 



TABLE 5 

Combined-Spectrum Cross-Correlation 



Major Minor 



a 




height 


velocity 2 


height 


velocity 


0.00- 


1.50" 


0.68 


+0.01 ±0.08 


0.67 


-0.01 ±0.08 


0.50- 


1.50" 


0.26 


-0.29 ±0.33 


0.28 


-0.42 ±0.28 


0.00- 


0.50" 


0.70 


-0.09 ±0.08 


0.69 


-0.03 ± 0.08 


0.50- 


1.00" 


0.23 


-0.17±0.41 


0.27 


-0.56 ±0.32 


1.00- 


1.50" 


0.16 


-0.17 ±0.27 


0.16 


-0.39±0.31 



"Velocities in km s 



TABLE 6 
Rotation Velocities" 



a 




Major 


Minor 


0.00- 


0.50" 


-0.03 ± 0.04 


-0.03 ±0.05 


0.50- 


1.00" 


-0.04 ±0.09 


-0.05 ±0.12 


1.00- 


1.50" 


-0.12±0.11 


-0.12 ±0.14 



"Velocities in km s 



which are most consistent with all correlations. In other 
words, cross-correlating the spectrum extracted from r, and 
-r, will give a measure of 2v,. However, one can also corre- 
late spectra from different radii, r, and rj, to obtain the relative 
velocity difference Av;;. 

Using all six narrow extractions yields 15 correlation per- 
mutations. Thus we have 15 measurements which depend on 
only three quantities, assuming that we should find equal but 
opposite velocities at the same radius on either side of the 
cluster. We solve for the most probable velocities by minimiz- 
ing the sum of the squares of the deviation of each correlation. 

Thus instead of basing a rotation velocity on only one mea- 
surement of the cross-correlation of two spectra at equal dis- 
tances from the cluster center, we measure 15 correlations be- 
tween all 6 spectral extractions and find the three velocities 
which are in best agreement with all measurements. These 
"best fit" rotation velocities are given as a function of semi- 
major axis distance in Table [6] We find no statistically sig- 
nificant rotation along the line of sight at any distance along 
either axis of WLM- 1 . 

5.1. Rotation Errors 

In order to estimate the robustness of our results we have 
performed various cross-correlation simulations. All of these 
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FIG. 7. — Measured cluster rotation along the major (top) and minor (bot- 
tom) axes. The points are taken from the last three lines of Tablel5l(divided 
by two to convert to rotation velocity). The data shown by the error bars are 
the results of simultaneously solving for the best rotation velocities using all 
15 narrow extraction cross-correlation permutations (Tabled. 



tests start out with a base spectrum which is the average of the 
major- and minor-axis integrated spectra of WLM-1. 

To test the dependence of the cross-correlation velocity on 
the signal-to-noise ratio of the spectra we create pairs of spec- 
tra at various signal levels by adding random noise to two 
copies of the base spectrum, treating them as opposite sides 
of the cluster. We then remove the telluric absorption and 
shift one of the spectra by 5 km s" 1 . The pair of spectra 
are then cross-correlated as described in Section [5] to obtain 
their relative velocity. The results of these tests are illustrated 
in Figure [8] along with the height and FWHM of the cross- 
correlation peak. This figure shows no systematic trends with 
signal-to-noise ratio, and verifies that the velocity errors re- 
turned by IRAF are indeed valid. From this analysis we es- 
timate that the errors associated with the measurement of the 
cross-correlation FWHM are approximately three times the 
reported velocity errors, and have plotted them as such in Fig- 
ureEl 

To estimate the errors of the x 2 -minimization technique 
used to determine the "best" rotation velocity we ran 50 sim- 
ulations consisting of six copies of the integrated spectrum of 
WLM-1 to mimic the six 0.5" extractions. Each of the six 
spectra is shifted such that the simulated cluster has a mono- 
tonically increasing rotation velocity with radius, and equal 
but opposite rotation velocities on opposite sides of the clus- 
ter. Specifically, the spectra were shifted by -3, -2, -1, +1, 
+2, and +3 km s" 1 each. Noise was added to each spectrum 
so that the cross-correlation height measured in the simula- 
tion (Fig. [HJ matched that measured in the real data (columns 
2 and 4 in Table lU. 

The velocity error distributions (recovered - input) result- 
ing from the 50 simulations have widths of 0.02, 0.05, and 
0.06 km s" 1 for the central, middle, and outer velocities re- 
spectively. For each simulation we also record the average er- 
ror of each measurement which results from making all cross- 



correlation measurements consistent (the square of this quan- 
tity is actually what is minimized when solving for the best 
set of rotation velocities). The mean of this error over all sim- 
ulations is 0.031 km s . However, in the real data, the error 
is 0.056 km s" 1 along the major axis and 0.074 km s" 1 along 
the minor axis. This discrepancy indicates that the simulation 
was too idealistic, and we have therefore scaled the estimated 
velocity errors by 1.8 for the major axis, and 2.4 for the mi- 
nor axis, yielding our formal error estimates of 0.04, 0.09, 
and 0.1 1 km s" 1 for the central, middle, and outer extractions 
along the major axis, and 0.05, 0.12, and 0.14 km s" 1 along 
the minor axis. 




50 100 
Additional Read Noise 



150 



FIG. 8. — Results of cross-correlation simulations using the combined 
major- and minor- axis integrated spectra with varying amounts of additional 
random read noise. The spectra were offset by 5 km s , and here we illustrate 
the difference in the cross-correlation velocity (top), and the width (middle) 
and height (bottom) of the cross-correlation peak. The velocity errors are 
those returned by FXCOR, and the FWHM errors are three times the velocity 



6. CLUSTER VELOCITY DISPERSION 

The velocity dispersion can be derived from the width of 
the cross-correlation peak, which is the sum of the line widths 
of the two spectra being correlated added in quadrature. In the 
case of the globular cluster WLM-1, the line width is a combi- 
nation of the intrinsic stellar line width, the instrumental pro- 
file, and the velocity dispersion of the cluster. Unfortunately, 
no template star was measured at the time of our observations. 

In order to estimate the velocity dispersion of WLM-1, we 
used two high-resolution stellar sp ectra obtained from the 
UVES Paranal Observatory Project (Baenul o et al.l l2003) as 
templates. We chose HD 30562 and HD 45067 (both F8V 
stars) because of their similarities in line strengths to the 
WLM-1 spectrum. It is important to note that the template 
spectra were obtained using the 0.5" slit and no binning, while 
the WLM- 1 spectra used 2x2 binning. We therefore rebinned 
the template stars to the same resolution as our observations. 

In order to increase the signal-to-noise ratio of the WLM-1 
spectra, we averaged together the spectra extracted from op- 
posite sides of the cluster (since we only want to compare 
the velocity dispersions along the major and minor axes). We 
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then analyzed the spectra measured on each chip separately 
to allow optimal Fourier filtering and rejection of undesirable 
lines. For example, we rejected the blue end of the blue chip 
which was contaminated by broad hydrogen lines and Call 
H&K lines, the region around H/3 on the Red-L chip, and Ha 
on the Red-U chip. 

We then cross-correlated each WLM-1 spectrum with each 
of the two template spectra. This gives a cross-correlation 
width for each template, which are then averaged. The result- 
ing cross-correlation widths derived for the integrated cluster 
spectra are 28.7 and 29.0 km s" 1 for the major and minor axes 
respectively. This width is the quadratic sum of the WLM-1 
and template line widths. In order to estimate the line widths 
of the templates, we cross-correlate the two templates with 
one another, and averaging the values from each chip, we find 
a width of 16.4 km s" 1 . 

If we make the (generous) assumption that the WLM-1 and 
template spectra have the same instrumental width and intrin- 
sic stellar line widths, it is now possible to subtract these com- 
ponents to find the stellar velocity dispersion. Subtracting in 
quadrature yields velocity dispersions of 23.5 and 23.9 km s" 1 
along the major and minor axes respectively. Following the 
same procedure for each of the various spectral extractions of 
WLM-1 yields the velocity dispersion profile shown in Fig- 
ure |9] where the err ors a re three times the velocity errors, as 
estimated in Section lSTl Here the extraction semi-major axis 
distances are shown using flux-weighted centers. 

For comparison we have overplotted the velocity dispersion 
profiles predicted by the isotropic Maxwellian distribution of 
the King model. These curves sho w the velocity dispersions 
computed according to eq. (31) in iKind (Q966) with a con- 
centration c = 2.00, 2.25, and 2.50, scaled to m atch the ob- 
served central velocity dispersion (but see Sec. 16. 1> . where 
the c = 2.5 profile is the flattest and, of these three curves, 
the best fit to the data. The discrepancy between the model 
and the measured data points (given that the light profile indi- 
cates c = 2.1 ±0.3) suggests that either the mass distribution is 
more extended than predicted by a King profile or that there 
is a systematic problem in our measurement of the velocity 
dispersion. 

In order to verify that any observed differences (or lack 
thereof) in line widths are inherent in the source spectra and 
not due to changes in the instrumental/atmospheric contribu- 
tion, we have also compared the widths of unresolved telluric 
absorption lines in each using the atmospheric transmission 
spectrum described in Section l4~2l as the template. We first 
excise a region of each WLM-1 spectrum in the range 6270- 
6320 A which is dominated by absorption lines due to atmo- 
spheric O2. We remove the continuum and broad features 
by fitting with a high-order cubic spline, and estimate the in- 
strumental width of each via cross-correlation with the atmo- 
spheric spectrum. The widths of the telluric lines present in 
the major- and minor-axis spectra differ by less than 0.05 km 
s" 1 , confirming the stability of the conditions during the ob- 
servations. 

Our absolute measurement of the velocity dispersion pro- 
files is uncertain due to the low signal-to-noise of our spec- 
tra, and the imperfect match with the templates which were 
obtained at a different time using a different instrument con- 
figuration. However, the relative velocity dispersion between 
the major and minor axes should be much more reliable. The 
bottom of Figure [9] shows the difference between the veloc- 
ity dispersions measured along the major and minor axes. 



All spectral extraction regions are consistent with zero differ- 
ence, and comparing the integrated cluster spectra measured 
along the major and minor axes yields a difference of only 
-0.4 ± 2.6 km s" 1 (major-minor). 




0.5 1 1.5 

semi— major axis (arcseconds) 



FIG. 9. — The velocity dispersion estimated via cross-correlation along the 
major-axis (filled circles) and minor-axes (open circles) of WLM-1 (top). The 
curves show the velocity dispersion profiles predicted by a King model with 
c = 2.00, 2.25, and 2.50. The bottom panel shows the difference in velocity 
dispersion between the major and minor axes. 



6.1. Central Velocity Dispersion Estimate 

We can obtain an estimate of the central velocity dispersion 
of WLM-1 as a reality check of our measurement of the ve- 
locity dispersion profile. The simplest way to obtain a rough 
estimate of the velocity dispersion is to utilize the correlation 
between the central velocity dispers ion and the cluste r abso- 
lute magnitude which was found by Diorgo vski et alJ (2003) 
to be followed by globular clusters in several different galax- 
ies. On the basis of Figure 3 (top panel) in Djorgovski et al., 
we find 

logcr = -1.01 -0.229M y , (1) 

where a is in km s" 1 . We use the apparent magnitud e 
of the cluster from Table I in San dage & Carlsonl \ 19851) . 
V = 16.06 mag, and the extin ction and distance moduli esti- 
mated by Hodg e et alJ J1999I) . Ay = 0.07 mag and (m-M) = 
24.73 mag, respectively, to find My = -8.74 mag for WLM- 
1 . The above relation then gives a central velocity dispersion 
cr~ 9.810ns- 1 . 

A perhaps more accurate estimate for the central velocity 
dispersion of the cluster can be obtained from the "funda- 
me ntal plane" correlation for globular clusters, as provided 
by Diorgovski (1995). In particular, Djorgovski has shown 
that (Galactic) globulars follow tight relations between the 
central velocity dispersion er, the central surface brightness 
/j,v{0), and the core radius r c , on the one hand; and between 
er, the half-light radius r/,, and mean surface brightness within 
r/,, (nv)h- Based on his equations (3) and (6), we find: 
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log a = 4. 173(±0. 175) +0.45 log r c - 0.204(±0.008) MO), 

(2) 

and 



log cr = 4.829(±0.237) + 0.70 log r h - 0.244(±0.012) {fj, v ) h , 

(3) 

where r c and r/, are in parsecs, and a in km s 1 . 

Usin g equ ation (|2) with the structural parameters derived in 
Section l2~2l core radius r c = 0.56 pc and central V-band sur- 
face brightness nv(0) = 14.83 mag/arcsec 2 , we estimate the 
velocity dispersion to be lO.S^ 5 km s" 1 . Equation (|3j with 
the half-light radius r/, = 3.21 pc and the mean surface bright- 
ness inside the half-light radius (/ly)/, = 17.31 mag/arcsec 2 
gives a central velocity dispersion of 9.1+5 6 g 3 km s . These 
estimates of the central velocity dispersion are lower than our 
measurements in Figure |9] but the uncertainties do not ex- 
clude our value. The larger velocity dispersion implies that 
the mass of WLM-1 lies on the high side of the relation de- 
rived for Galactic clusters. 

7. CLUSTER VELOCITY 

Extracting an integrated spectrum from the entire cluster 
using a 3" (±1.5") aperture yields a relatively high signal- 
to-noise spectrum relative to the partial extractions used to 
measure cluster rotation. These integrated spectra have many 
strong lines from which we can measure the net radial velocity 
of the cluster. In order to make this measurement, we first 
calculated and corrected for terrestrial motion, based on the 
observation times of 2003/10/03 at 3:23 and 4:23 UT for the 
minor and maj or-axis respectively. W e obtained an atomic 
line list from iReader & Corliss! d!980h . and used the IRAF 
task RVIDLINES to calculate the radial velocity. Centroiding 
~ 70 lines in each of the major and minor axis spectra yields 
heliocentric radial velocities of - 1 05 .77 ± 0.50 and - 1 05 .84 ± 
0.50 km s" 1 respectively. Averaging the results we estimate 
that the heliocentric radial velocity of WLM-1 is -105.8 ± 
0.4km s" 1 . The radial velocity of the WLM galaxy is -1 16 ± 
2 km s" 1 (Huch ra et alJ fl993). thus the cluster has a velocity 
of 10 km s" 1 with respect to its host. 

8. VELOCITY ANISOTROPY 

It appears that the ellipticity of WLM- 1 cannot be due to ro- 
tation, so perhaps WLM-1 is (in a sense) more like elliptical 
galaxies, where anisotropies in the velocity dispersion play 
an important role in the observed morphology. To address 
this hypot hesis we employ the tensor virial theorem, as de- 
scribed in iBinnev & Tremainel dl987l) . to obtain information 
about the relation between the ratio of rotation velocity and 
velocity dispersion v rot /a, on the one hand, and ellipticity, on 
the other. 

We start by assuming that a globular cluster is an oblate 
spheroid, rotating around the z axis, and is seen edge-on. We 
assume, moreover, that the system's isodensity surfaces are 
concentric ellipsoids. Defining v 2 ot0 as the mass-weighted 

mean-square rotation speed, and <r ( 2 as the mass-weighted 
mean-square random velocity along the line of sight, we ar- 
rive at the following relation: 
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FIG. 10. — The predicted ratio between mass-weighted rotation velocity 
and mass-weighted velocity dispersion is plotted as a function of the elliptic- 
ity for an oblate spheroid. The curves correspond to values of the velocity 
anisotropy factor S, with an interval of 0. 1 in 8 from the upper to the lower 
curve. The dashed line indicates the measured ellipticity value for WLM- 1 . 



where S < 1 is a parameter that measures the anisotropy of 
the velocity dispersion tensor and e is the eccentricity, which 
is related to the ellipticity e as follows: 



Vl-d-e) 2 - 



(5) 



Using these equations, we can build a family of curves giv- 
ing v rot o /(To as a function of ellipticity for different values of 
6. The result is shown in Figure ^3 Binney & Tremaine 
(1987) argue that the representative points in the (v rot /er,e) 
plane of systems that are flattened by rotation lie close to the 
curve (5 = 0, independent of the inclination angle i. 

The mean ellipticity of WLM-1, as measured in Section l2~Tl 
is e = 0.17, which is indicated in Figure^] by the vertical 
dashed line. As the low mass of the WLM galaxy effectively 
rules out the possibility that WLM-1 owes its ellipticity to 
tidal effects or gravothermal shocks, this figure very nicely 
illustrates the kinematical combinations available to explain 
the ellipticity of WLM-1. 

At one extreme, the ellipticity of WLM-1 could be pro- 
duced entirely by rotation. In this scenario there is no ve- 
locity anisotropy (5 = 0), and reading from the plot, we see 
that v ro t/er ~ 0.57. Taking a conservative estimate of 10 km 
s" 1 for the central velocity dispersion of the cluster (Sec. |6]», 
the required rotation would be 5.7 km s" 1 . However, as we 
find no measurable rotation (v lot < 0.1 km s" 1 ), we exclude 
this possibility. 

At the other extreme, the ellipticity of WLM-1 is due en- 
tirely to velocity anisotropy. The quantity (1-8), defined as 
the ratio between the diagonal terms of the velocity dispersion 
tensor along the rotation ("z") axis and along the elongation 
("jc") axis, should scale as 



(4) 



(6) 



From Figure^] one sees that if there is zero rotation (v rol /<T = 
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0), it is possible to produce the observed ellipticity with a ve- 
locity dispersion parameter of <5 = 0. 14. Therefore, WLM- 
1 's flattening could be accounted for if the velocity dispersion 
along the rotation axis were smaller than along the elongation 
axis by 

a z xs 0.93 a z . (7) 

Again taking the nominal velocity dispersion to be ~ 10 km 
s , this velocity anisotropy would only produce a difference 
of 0.7 km s" 1 in velocity dispersion between the major and 
minor axes. Such a difference would be too small to be de- 
tected on the basis of our data. Therefore, while we cannot 
prove that velocity anisotropy is responsible for the observed 
shape of WLM- 1 , it certainly appears to be the only option 
that is consistent with our observations. 

9. CONCLUSIONS 

We have analyzed VLT UVES long-slit spectra and HST 
WFPC2 images of the extragalactic globular cluster WLM-1 
in an attempt to understand the cause of its elliptical morphol- 
ogy. Unlike the situation for Galactic globular clusters whose 
detailed shapes may be influenced by their massive host 
galaxy, WLM is a low-mass dwarf irregular galaxy, eliminat- 
ing the possibilities that the non-sphericity of WLM-1 could 
be caused by either tidal stresses or gravothermal shocks. This 
leaves two explanations for the flattening of this cluster: in- 
ternal rotation and/or velocity dispersion anisotropy. 

Cluster rotation seemed to be a prime candidate for the ori- 
gin of WLM-l's ellipticity since two-body relaxation, which 
erases any primordial anisotropy, is relatively efficient in star 
clusters. In addition, rotation has been measured for several 
Galactic globular clusters using bot h proper motion and radial 
velocity techniques [the review by M evlan & Heggid i 19971) 
lists the rotation velocities of 11 clusters]. The most notable 
case is that of uj Cen, which exhibits correlated rotation and 
flattening profiles, strongly suggesting that its ellipticity is due 
to rotation. 

We therefore obtained high-resolution long-slit echelle 



spectra along the major and minor axes of the cluster. We used 
cross-correlation to look for velocity shifts on opposite sides 
of the cluster, dividing the cluster into annuli in an attempt to 
obtain a rotation velocity profile. However, we found no ev- 
idence for rotation (v 10t < 0.1 km s" 1 ) along either the major 
or minor axes of the cluster. As the cluster would require a 
rotation of > 6 km s" 1 to sustain its shape, we conclude that 
rotation is definitely not the main driver of WLM- 1 's elliptic- 
ity. 

The only remaining plausible source of flattening in the 
WLM cluster is velocity dispersion anisotropy, the main rea- 
son elliptical galaxies are elliptical. We attempted to measure 
the stellar velocity dispersion in WLM-1 based on the width 
of the cross-correlation peak; however, due to our large errors, 
we could see no statistically significant difference between 
the velocity dispersions measured along the major and minor 
axes. Calculations reveal that the magnitude of the velocity 
anisotropy required to produce the observed ellipticity is rel- 
atively small: a 7% difference in velocity dispersion, or less 
than 1 km s -1 , which is smaller than our velocity dispersion 
measurement errors. 

In summary, we conclude that the ellipticity of the globular 
cluster WLM-1 is most likely due to anisotropy in the stellar 
velocity dispersion. 
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